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We study inhomogeneous states in the t-J model using an unrestricted Gutzwiller approximation. 
We find that pa x pa checkerboard order, where p is a doping dependent number, emerges from 
Fermi surface instabilities of both the staggered flux phase and the Fermi liquid state with realistic 
band parameters. In both cases, the checkerboard order develops at wave vectors (±27r/pa, 0), 
(0,±27r/pa) that are tied to the peaks of the wave-vector dependent susceptibility, and is of the 
Lomer-Rice-Scott type. The properties of such periodic, inhomogeneous states are discussed in 
connection to the checkerboard patterns observed by STM in underdoped cuprates. 



Spatially inhomogeneous, periodically modulated lo- 
cal density of states (LDOS) has been observed recently 
by STM in high-T c cuprates under a variety of condi- 
tions where the superconductivity is weakenediSiiik. 
The hope for a certain type of ordered state competing 
with superconductivity has been raised in connection to 
the pseudogap phenomena. The tunneling conductance 
map exhibits pa x pa checkerboard (CB) patterns corre- 
sponding to wave vectors (±2ir/p, 0), (0, ±27r/p), where 
p takes on values between 4 and 5 (hereafter, we set the 
lattice constant a = 1). In and around a vortex core in 
Bi 2 Sr2CaCu20 8+( 5 (BSCCO), p ~ 4.3i; p ~ 4 in opti- 
mally doped BSCCO in zero magnetic field 2 ; p ~ 4.7 in 
underdoped BSCCO in the pseudogap phase above T<£; 
p ~ 4 in lightly doped oxychlorides Ca2- a; Na ;E Cu02Cl2 
(Na-CCOC) at very low temperatures^; and in substan- 
tially underdoped BSCCO, p ~ 4.5 were observed in the 
dark regions of the conductance map^. The nature and 
the origin of the LDOS modulations, in particular the 
short coherence length and the role of dopants&i, are 
unclear at present. Several inhomogeneous electronic 
states have been proposed theoretically, including pair 
density waves&S; hole Wigner crystaUS; a Wigner crys- 
tal of hole pairs embedded in a d-wave resonating va- 
lence bond (RVB) state 11 ; and valence-bond solid with 
or without charge order—. How certain inhomogeneous 
electronic states with CB order arise from the micro- 
scopic t-J model of doped Mott insulators is the focus of 
this work. Recent density matrix renormalization group 
calculations performed on small t-J clusters found ap- 
proximate CB like patterns with strong one-dimensional 
stripe characters^. In this paper, we discuss a spa- 
tially unrestricted Gutzwiller approximation, which al- 
lows us to study large systems in the thermodynamic 
limit. We focus primarily on inhomogeneous solutions in 
the nonsuperconducting phase at moderate doping where 
the underlying homogeneous phase exhibits a Fermi sur- 
face (FS). Specifically, we consider the two situations 
displayed in Fig. 1 for the spectral intensity of the low 
energy single-particle excitations as measured by angle- 
resolved photoemission (ARPES), in the staggered flux 
phase (SFP)iA and the uniform short-range RVB state 15 
with realistic band parameters. The quasiparticle scat- 
tering in the SFP (Case A in Fig. la) is dominated by the 
wave vector q* connecting the tips of the Fermi pock- 
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FIG. 1: Spectral intensity near the Fermi level in the SFP (a) 
and the Fermi liquid phase with realistic band parameters (b) 
in the t-J model. Those in the corresponding checkerboard 
ordered states are shown in (c) and (d) where sections of the 
FS with high spectral intensity connected by q* are gapped. 



ets with high spectral intensity. In Fig. lb (Case B), 
the sections of the large FS around (±7r, 0) and (0, ±tt) 
are nested and scattering by the nesting vector q* is en- 
hanced. In both cases, the wave vector dependent sus- 
ceptibility exhibits sharp peaks at q* and the system is 
prone to superlattice instabilities, akin to the examples 
discussed by Lomer— , and Rice and Scott 17 for two- 
dimensional band structures with nesting and saddle- 
points respectively. We show that this indeed happens in 
both cases, leading to inhomogeneous CB ordered states 
that are lower in energy than the uniform state. The re- 
sulting spectral intensity maps are shown in Figs, lc and 
Id where FS sections connected by q* are truncated. In 
Case A, the CB order is a secondary instability of the SFP 
that already exhibits a pseudogap in the LDOS. Case B, 
however, is more significant. The CB order produces 
a pseudogap around (±7r, 0) and (0, ±71"), the antinodes 
of the d-wave pairing gap function, leaving behind only 
"Fermi Arcs", in agreement with ARPES experiments 18 . 

In the Gutzwiller approximation, the projection of 
double-occupation from the Hilbert space is partially ac- 
counted for by the statistical weight factors, g l and g J 
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multiplying the quantum coherent states connected by 
the hopping and the exchange processes^. The latter 
lead to renormalizations of the hopping and exchange 
parameters. The renormalized t-J model can be studied 
by decoupling the exchange term 
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Here /i is the chemical potential, Xij — ^}2a( c la c j^) anc ^ 
\j = ( c i] c ji — c ii c j]) are the the nearest neighbor bond 
and flux, and pairing fields, and tij is the hopping be- 
tween sites i and j. = t between nearest neigh- 
bors (NN); t' for next NN. In the last line of Eq. (J, 
Xi = 1— X^cr( c L c *cr) denotes the doped hole density at site 
z, and the averaged doping is given by S = (1/N S ) X{ 
on a lattice of N s sites. 

There are two new features in the unrestricted 
Gutzwiller approximation. First, the usual Gutzwiller 
factorsiS now depend on the local hole density, 
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In essence, the local doping concentration X{ is promoted 
to a variational parameteA Second, as an electron hops 
between sites, the renormalized bandwidth will change 
by 0(1/N S ). However, the kinetic energy of the occu- 
pied states changes by an amount of order unity, and 
this energy difference must be reflected in the equilib- 
rium condition by the local fugacity Si in the last term 
of Eq. QJ. The value of S{ is determined by minimizing 
the energy with respect to 
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In the rest of the paper, we focus on the normal state 
unless otherwise noted. 

Checkerboard order in the SFP We choose t = 3 J and 
set t' and further neighbor hopping to zero. In this case, 
the uniform normal state of Eq. is the SFP for small to 
moderate doping. The spectral intensity near the Fermi 
level is shown in Fig. la for 5 = 0.102. To locate the 
wave vectors at which Lomer- Rice- Scott instabilities may 
occur, we calculate the static susceptibilities 
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FIG. 2: (a) Susceptibility in the SFP in the density-bond 
channel along (q x ,0). The peak positions correspond to 
q* = (27r/4,0), q* = (2tt/5,0) for the two doping levels, (b) 
Density-bond susceptibility at fixed q* = (2tt/4, 0) (squares) 
and (2tt/5,0) (trangles) as a function of doping. Also shown 
is the doping dependence of inhomogeneity period p = 2ii/q* 
(circles). 



where k™ n (k,q) is the coupling vertex in the density, 
bond, and flux channels labeled by a, E% is the quasi- 
particle band dispersion^. We find that all susceptibil- 
ities exhibit peaks at q* = (g*, 0), (0, q*). The domi- 
nant density-bond (x-direction) susceptibility is plotted 
in Fig. 2a along (q x ,0). The q* can be identified by the 
location of the sharp peak and is doping-dependent and 
incommensurate in general. If an instability occurs at 
q*, the ground state will exhibit CB order with a pe- 
riodicity p = 2ir/q*, which is shown in Fig. 2b, as a 
function of doping. Note that the maximum periodicity 
is limited by the stability of the SFP. For our parame- 
ters, p varies continuously between 4 and 5 in the doping 
range of 10 — 12%. To search for a spatially unrestricted 
Gutzwiller solution that allows static order at q* in fi- 
nite systems, it is advantageous to study doping levels 
at which q* takes on commensurate values. In Fig. 2b, 
the susceptibility at fixed q* = (2tt/4,0) and (2tt/5,0) 
are shown as a function of doping. The peak structures 
allow the choice of doping to be made. We thus focused 
on an average doping S = 0.102 where q* is consistent 
with an instability toward p = 4 CB order. 

We next study spatially unrestricted solutions on 80 x 
80 systems with 4x4 and 8x8 supercells wherein the x\, 
Si, and Xij are allowed to have spatial variations. Their 
values are determined by the standard iterative solutions 
of the self-consistency equations. An additional four-fold 
symmetry may be enforced to reduce the number of vari- 
ational degrees of freedom for faster convergence. Re- 
markably, the self-consistent state with p = 4 CB order 
emerges with a significantly lower energy compared to 
the uniform SFP state. The real space CB patterns of 
local doping Xi, staggered plaquette flux (-lf^gX^ 
valence bond Xi an d integrated LDOS are shown in 
Fig. 3 (a-d). Their Fourier transforms show dominant 
peaks at q* = (±2tt/4, 0) , (0, ±2tt/4). The relative mod- 
ulations in the CB state and the energy comparison to 
the uniform state are summarized in Table I. Interest- 
ingly, the charge density variations are small (less than 
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2%) while the bond and the plaquette flux variations are 
more significant. The flux and bond density wave insta- 
bilities in the SFP of the t-J model were first revealed 
in the slave boson large-N theor^i. These results show 
that weak charge order in the SFP is a byproduct of the 
driving CB order of the valence bond and the plaquette 
flux. The accompanying weak charge ordering is in line 
with STM observations in Na-CCOGi 




FIG. 3: 2D maps of (a) local doping concentration, (b) stag- 
gered plaquette flux, (c) x-direction bond, and (d) integrated 
local density of states. 

It is important to note that the appearance of the local 
fugacity S{ in Eq. Q is essential for the emergence of the 
CB state with lower energy than its uniform counterpart. 
Indeed, without Si, a much more inhomogeneous state 
can emerge as a self-consistent solution^, but its energy 
is significantly higher than the uniform state (see Table 



TABLE I: Modulation amplitudes and the energy of the 4x4 
and the 5x5 checkerboard states emerging from the SFP and 
the Fermi liquid state. 
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The projected spectral function in the Gutzwiller ap- 
proximation is given by A(i,j,ui) = —2Img t i jG(i,j,cu + 
i0 + ), where G is the real-space retarded Green's function. 
The Gutzwiller factor arises from projecting the matrix 
elements of the electron operator between ground and 
excited states. The spectral intensity at the Fermi level 
A(k, co ~ 0) maps out the FS shown in Fig. lc. The 4x4 
checkerboard order opens and gaps out the tips of the 
Fermi pockets. The intensity is dominated by the inner 
branch of the remaining segments due to the anisotropy 
in the SFP coherence factors. Fine structures due to zone 



folding and higher order scattering effects become visible 
only after the intensity scale is reduced by four orders 
of magnitude. In Fig. 4a, the band dispersion E(k) ex- 
tracted from the quasiparticle peaks in A(k,w) is shown. 
It retains the general structure of the SFP, exhibiting a 
pseudogap in the LDOS p(r,u) = ^ k e ikr A(k,cu) shown 
in Fig. 4b. The particle-hole asymmetry is much en- 
hanced at these parameters. The checkerboard order in- 
duces a small energy gap at Ep in the dispersion along 
the M-X direction, which is reflected in the additional 
dip in the LDOS at the Fermi level. It is worth point- 
ing out that the suppression of the LDOS is particle-hole 
asymmetric and predominately resides on the occupied 
side. 
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FIG. 4: (a) Quasiparticle dispersion in the checkerboard SFP. 
Inset shows the gap near the FS. (b) LDOS site (1, 1). Inset 
shows the LDOS near Fermi energy. 



Checkerboard order with realistic band dispersion Next 
we examine whether the CB state emerges in the t-J 
model with further neighbor hopping capable of produc- 
ing the dispersion and the FS observed by ARPES 22 . 
The hopping integrals in Eq. are chosen to be 
6.0,-2.0,0.71,0.7,-0.41,0.07 in units of J from the 
nearest up to the sixth neighbors respectively. The cor- 
responding FS in the uniform Fermi liquid was shown in 
Fig. lb at doping 5 = 0.11. It reveals the partially nested 
segments near (±7r, 0) and (0, ±7r) with nesting vectors 
q* ~ (±27r/5,0), (0, ±27r/5). The set of parameters is 
chosen such that the nesting properties in the realistic 
band dispersion are most pronounced at low dopings. 
Similar to the SFP, the calculated static, wave vector 
dependent susceptibilities peak at these nesting vectors 
and are most pronounced in the density-bond channel. 
Fig. 5a displays the (logarithmic) divergence as q — > q*, 
indicating possible instabilities toward spatially inhomo- 
geneous states. On 80 x 80 systems with 10 x 10 su- 
per cells, we find that the spatially unrestricted solution 
converges to the CB state with p = 5, consistent with 
the nesting vector q*. The 5x5 pattern is shown in 
Fig. 5b for the doping concentration. The amplitudes of 
the variations in doping and bond (Table I) are signifi- 
cantly larger than in the SFP owing to the more singular 
behavior of the susceptibility. The quasiparticle disper- 
sion obtained from the calculated A(k, uo) near the Fermi 
level is shown in Fig. 5c. The CB order opens up a siz- 
able gap near (±7r, 0), (0, ±tt) and truncates the nested 
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FIG. 5: (a) Density-bond susceptibility along (q Xj 0). (b)2D 
map of the local doping showing 5a x 5a order, (c) Quasi- 
particle dispersion, (d) Tunneling LDOS. The circle marks 
the gap opening just below the Fermi level. Inset shows the 
spectra near Fermi level at all different sites. 



segments of the original FS as shown in Fig. Id. 

The gap near the antinodes and the residual Fermi arcs 
around the nodes are reminiscent of the pseudogap phe- 
nomenology observed by ARPES 23 . It is important to 
point out the differences between the CB order induced 
pseudogap and one that is due to a d-wave pairing gap. 
A finite zero-bias conductance must remain due to the 
finite "volume" of the residual Fermi arcs, whereas a d- 
wave paring gap leads to a linearly vanishing LDOS. In 
Fig. 5d, the calculated p(r, uo) in the p = 5 CB state is 
shown. The gap opening pushes the van Hove peak down- 
ward, resulting in a substantial lowering of the kinetic 
energy. However, the gap opens predominantly on the 
occupied side leaving the spectrum particle-hole asym- 
metric even at low energies. We find that this feature, 
also observed in the inhomogeneous SFP in Fig. 4b, is 



rather generic of the CB ordered states associated with 
the partial gapping of the FS. This is in contrast to the 
particle-hole symmetric LDOS in a pairing-induced pseu- 
dogap state and the conventional charge density wave 
state with a fully gapped FS. A direct comparison to the 
STM spectra in oxychlorides 4 is difficult since the cor- 
relation length of the observed CB pattern is very short 
(20a), although particle- hole asymmetry in the tunneling 
spectra is apparent. 

In conclusion, we have shown that pa x pa ordered 
states emerge in the t-J model from both the SFP and the 
Fermi liquid phase with realistic band parameters for the 
cuprates. The CB states we found originate from a gener- 
ally incommensurate partial FS instability of the Lomer- 
Rice- Scott variety, and are different from the bond- 
ordered states 12 inherent of the spin-Pierls or dimerized 
spin liquids at half-filling. As such, the inhomogeneous 
states found in the unrestricted Gutzwiller approxima- 
tion are expected to be robust against Gutzwiller projec- 
tion of double occupation. We have tested that the CB 
state is stable in the presence of a weak pairing order pa- 
rameter and long-range Coulomb interaction. However, 
in the parameter regime investigated, the uniform d-wave 
superconducting state always has a lower energy. Thus, 
the CB ordered state cannot coexist with d-wave pairing 
well inside the superconducting phase where the FS has 
been gapped out except for point zeros (nodes). Due to 
such competitions, the CB order discussed here may arise 
in the cuprates only in the pseudogap regime or close to 
the superconducting phase boundary, and may manifest 
in the form of fluctuating order pinned by dopant dis- 
order. Finally, the NN Coulomb interaction V in the 
t-J-V model, which is known to reduce the d-wave pair- 
ing strength while enhancing that of the valence bond, 24 
may lead to the coexistence of CB order and d-wave su- 
perconductivity. 
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